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Abstract 


The  thermally  grown  InP  oxide  as  etched  by  an  aqueous  dilute 
HF  solution  has  been  studied  by  ellipsometric  techniques.  The  ex- 
situ  measurement  reveals  a  two-layer  structure  for  the  oxide  grown 
at  440°C.  The  refractive  indices  for  both  oxide  layers  have  been 
determined  using  a  two-layer  optical  model.  The  etching  process 
has  also  been  monitored  ellipsometrically  in  the  real  etching 
environment,  in-situ.  A  fused  silica  cell,  which  enables  the 
windows  to  be  aligned  properly,  has  been  specifically  designed  for 
the  in-situ  solution  measurement.  A  liquid  layer  at  the  solution- 
oxide  interface  has  been  identified,  and  the  layer  is  shown  to 
contain  P  and  In  species  resulting  from  the  etching  reactions.  A 
theory  based  on  the  Lorentz-Lorenz  relation  results  in  a  reasonable 
qualitative  description  of  the  liquid  layer.  During  the  etching  of 
the  oxide  the  liquid  layer  shrinks  at  a  linear  rate,  and  after 
removal  of  the  outer  oxide  layer  the  liquid  layer  forms  a  dense 
electric  double  layer. 
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Introduction 


The  surface  oxide  layer  on  semiconductor  surfaces  continues  to 
attract  considerable  research  effort.  The  recent  applications  in 
opto-electronics  and  microwave  devices  have  drawn  increasing 
attention  towards  InP.  The  oxide  thermally  grown  on  InP  has  been 
widely  studied  in  an  attempt  to  use  it  as  an  insulating  material 
for  the  InP  surface (1-7) ,  as  has  been  successfully  accomplished  for 
Si02  films  on  silicon.  The  chemical  composition  of  the  oxide  has 
been  found(5,7)  to  be  both  complicated  and  process  dependent. 
Although  some  questions  still  remain,  it  is  generally  accepted  that 
the  oxide  grown  at  low  temperature  (T  <  600°C)  is  composed  of  two 
layers  of  differing  chemical  composition:  an  In-rich  outer  layer  of 
a  mixture  of  ln203  and  InP04  and  an  inner  layer  of  predominantly 
InP04(8,9).  Recently,  spectroscopic  ellipsometry  has  been  applied 
(9,10)  to  examine  the  dielectric  response  of  the  InP  thermal  oxide, 
and  the  dielectric  function  in  the  energy  range  2-5  eV  for  the 
oxide  thermally  grown  at  several  temperatures  has  been  reported. 

Traditionally,  the  electric  double  layer  formed  at  the  solid 
aqueous  solution  interface  has  been  studied  using  capacitance 
voltage  measurements,  potentiometric  titrations,  and  microelectro¬ 
phoresis  methods (11) .  There  are  several  investigations  of  the 
specific  ion  adsorption  and  the  electric  double  layer  using 
ellipsometric  techniques (12, 13) .  Several  studies  have  focused  on 
crystalline  and  amorphous  Si02  immersed  in  an  aqueous  solution 
(14,15).  Recently,  Gould  and  Irene(16)  have  used  in-situ 
ellipsometry  to  study  the  Si02  surface  behavior  dynamically  by 
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monitoring  changes  in  thicknesses  and  refractive  indices  of  the 
oxide  caused  by  a  reactive  solution,  i.e.  dilute  aqueous  HF, 
buffered  HF,  and  aqueous  NH4OH/NH4F,  and  Law  (17)  has  reported 
observations  of  orientational  ordering  of  water  and  organic 
solvents  on  pyrex  surfaces  by  in-situ  ellipsometry. 

In  this  paper,  we  report  an  ex-situ  ellipsometric  study  of  the 
etch  of  an  InP  oxide  thermally  grown  at  440°C  by  a  dilute  aqueous 
HF  solution.  We  show  clear  evidence  for  a  two-layer  structure  for 
the  oxide,  and  values  for  refractive  indices  for  the  outer  and 
inner  layer  of  the  oxide  at  the  wavelength  of  632. 8nm  have  been 
determined  using  a  two-film  optical  model.  We  also  describe  an  in- 
situ  observation  of  the  same  etching  system.  A  liquid  layer,  which 
has  a  different  optical  response  than  the  bulk  solution,  has  been 
found  to  be  present  at  the  oxide-solution  interface  during  the 
oxide  etching,  and  a  microscopic  model  based  on  the  Lorentz-Lorenz 
relation  has  been  used  to  describe  the  liquid  layer.  In  addition, 
we  report  the  observation  of  an  electric  double  layer  at  the 
solution-oxide  interface  which  forms  when  the  outer  oxide  layer  is 
removed. 

Experimental  Procedures 

N-type  undoped  (1  0  0)  InP  (n<l. 0xl016cm‘3)  wafers  were  cut  into 
lxl  cm2  samples.  The  samples  were  ultrasonically  degreased  in 
boiling  trichloroethylene  for  10  minutes,  followed  by  a  sequence  of 
acetone,  deionized  water (d.  i.  water)  rinse.  After  being  dipped  in 
a  concentrated  aqueous  HF  solution  for  15  seconds,  rinsed  in  d.  i. 
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water,  and  dried  in  N2  gas,  the  samples  were  loaded  onto  a  fused 
silica  boat.  Then,  the  boat  was  placed  in  the  end  cap  of  the 
oxidation  furnace  for  further  drying,  pushed  into  the  hot  zone  for 
a  15  minute  pre-oxidation  annealing  at  440°C  in  a  N2  atmosphere, 
oxidized  in  1  atm  02  at  440°C  for  4  hours,  and  annealed  in  N2  at  the 
same  temperature  for  15  minutes. 

A  commercially  available  manual  ellipsometer  automated  in  our 
laboratory  and  shown  in  Fig.  1  consists  of  a  HeNe  laser,  the  light 
source,  a  quartz  Rochon  polarizer,  a  three-dimensionally  adjustable 
sample  stage,  an  autocollimating  and  alignment  telescope (ACAT) ,  a 
rotating  analyzer,  and  a  signal  detection  system.  The  ACAT  is  used 
for  aligning  the  ellipsometer  itself  and  samples  to  be  measured. 
A  Gian  Tayor  polarizing  prism  that  is  mounted  on  a  hollow  motor- 
driven  shaft  serves  as  the  rotating  analyzer.  As  the  shaft 
rotates,  an  optical  encoder  generates  pulses  which  trigger  the 
analog  to  digital  conversions.  The  periodic  light  signal  through 
the  analyzer  is  detected  by  a  photomultiplier  tube(PMT) .  The 
output  PMT  current,  in  turn,  is  converted  to  voltage,  amplified  to 
0-10  volt  range,  and  converted  to  a  digital  signal  triggered  by  the 
encoder  pulses.  A  dedicated  personal  computer  controls  the 
measurements . 

A  fused  silica  cell,  as  shown  in  Fig.  2,  has  been  designed  for 
the  in-situ  etching  experiments.  Two  optically  flat  fused  silica 
plates,  serving  as  the  entrance  and  exit  windows,  are  connected  to 
the  main  part  of  the  cell  through  ball-socket  joints,  which  enable 
the  windows  to  be  adjusted  and  then  securely  set  in  place.  The 
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windows  have  to  be  orthogonal  to  the  incident  light  beam  for  the 
solution  measurements,  otherwise  the  media  with  different  indices 
of  refraction  on  each  side  of  the  windows  will  alter  the  direction 
of  the  light  beam,  hence  the  angle  of  incidence,  when  the  light 
beam  passes  through  the  windows.  The  solution  inlet  and  outlet 
permit  the  changing  of  a  solution  without  exposing  the  sample  to 
air. 

With  a  conventional  analyzer  in  the  position  of  the  rotating 
analyzer,  the  ellipsometer  and  its  optical  components  were  aligned 
and  calibrated  as  previously  described (18) .  Then,  the  conventional 
analyzer  was  removed  from  the  bench,  and  the  rotating  analyzer 
installed.  The  bench  was  set  in  a  straight  through (180°)  position 
and  a  series  of  measurements  were  carried  out,  and  the  offset 
correction  for  the  rotating  analyzer  was  determined  from  these 
measurements.  Then,  the  incident  angle  was  set  to  70.00°  for  all 
the  ellipsometry  measurements. 

The  etchant  used  in  both  ex-situ  and  in-situ  measurements  was 
1:500  HF(49%) :H20  solution.  For  the  ex-situ  measurements,  the 
samples  were  dipped  in  the  etchant  for  a  certain  time,  rinsed  in  d. 
i.  water,  and  blown  dry  by  N2  gas.  Prior  to  an  ellipsometric 
measurement,  the  sample  was  placed  into  a  slot  on  the  sample  stage 
so  that  the  same  surface  spot  could  be  probed  for  the  successive 
measurements,  and  the  sample  alignment  was  checked  and  adjusted  if 
necessary.  This  process  was  repeated,  with  the  etching  times  and 
ellipsometric  parameters  recorded. 

The  in-situ  measurement  first  requires  the  alignment  of  the 
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solution  cell.  An  sample  was  affixed  with  wax  to  a  fused  silica 
flat  inside  the  cell,  which  was  mounted  on  the  adjustable  stage. 
With  both  windows  absent,  the  sample  was  aligned (19) .  Then,  with 
the  windows  in  place,  the  cell  was  filled  with  d.  i.  water,  and 
the  entrance  window  was  aligned  by  adjusting  the  reflected  beam 
from  the  window  to  coincide  with  the  incident  beam,  and  then  the 
exit  window  was  aligned  independently.  The  first  measurements 
after  alignment  were  carried  out  in  d.  i.  water  to  insure  correct 
alignment,  and  then  the  dilute  HF  solution  was  rapidly  introduced 
into  the  cell  through  the  inlet  while  the  water  drained  through  the 
outlet.  During  the  entire  experiment,  fresh  HF  solution  was 
circulated  into  the  cell  at  a  rate  of  7-9  ml/min.  At  the  end  of 
the  experiment,  the  cell  was  drained,  thoroughly  washed  with  d.  i. 
water,  dried  by  N2  gas,  and  several  ellipsometric  measurements  were 
carried  out  in  the  air  ambient. 

Results  and  Discussion 

A  four  phase  model  consisting  of  ideal  optically  isotropic  and 
homogeneous  phases:  semi-infinite  ambient  -  film  1  -  film  2  - 
substrate,  has  been  used  in  the  ellipsometric  data  analysis  for 
both  ex-situ  and  in-situ  experiments (20) .  For  this  optical  system, 
there  are  three  interfaces  for  reflection  and  two  homogeneous  and 
isotropic  transmission  layers  involved.  The  reflection  and 
transmission  properties  for  each  interface  and  each  layer  can  be 
expressed  individually  by  a  characteristic  2x2  Jones  matrix,  and 
the  overall  properties  of  the  multilayer  structure  can  be  obtained 
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from  the  product  of  these  matrices (20) .  The  relationship  between 
the  ellipsometric  parameters  tfandA,  where  is  the  change  in  the 
amplitude  ratio,  and  A  the  change  in  the  phase  difference  between 
the  p  and  s  components  of  a  monochromatic  polarized  plane  wave  upon 
reflection  from  a  surface,  and  the  optical  properties  of  the  system 
is  given  by: 

tan  ¥exp(jA)  =  p( N^N^L^N^L^N^A)  (1) 

where  N0,  Nv  Nz,  and  Ns  are  the  complex  indices  of  the  ambient,  film 
1,  film  2,  and  the  substrate,  respectively.  L1  and  are  the 
thicknesses  for  the  films.  <j>  is  the  angle  between  the  incident 
beam  and  the  normal  to  the  sample  surface,  and  A  is  the  radiation 
wavelength.  With  the  N0,  <f>  and  A  as  experimentally  controlled 
variables,  and  with  Ns  fixed,  with  either  N,  and  L1  or  N2  and 
known,  the  other  pair  can  be  determined  from  the  measured  'PandA. 

Later  in  the  "In-Situ  Ellipsometry  Studies"  part,  we  model  a 
liquid  layer  as  a  medium  with  various  ions  and/or  molecules  near 
the  oxide  surface.  For  this  purpose  we  assume  a  Lorentz-Lorenz 
relation  which  treats  a  film  as  composed  of  polarizable  points. 
The  measured  refractive  index  can  be  related  to  the  microscopic 
properties  of  molecules  and  ions  which  constitute  the  optical 
system  by  the  Lorentz  Lorenz  equation (21) .  For  a  single  component 
system,  the  equation  can  be  expressed  as: 


4  n  a'  Na  ti  _  N1  -  1 
3  Mm  N*  ♦  2 


(2) 
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where  a’  is  the  polarizability  for  a  molecule  or  an  ion,  r?  is  the 
density,  and  M,,,  and  NA  are  molecular  weight  and  Avogadro's  constant, 
respectively.  With  the  assumption  that  the  microscopic  properties 
are  additive,  and  for  an  isotropic  and  homogeneous  phase  consisting 
of  m  non-interacting  components,  the  overall  optical  properties  can 
be  obtained  in  terms  of  the  optical  properties  of  the  individual 
components (22) : 


N1  -  1 
ff1  *  2 


ai 


(3) 


where  ^  and  a,  are  the  volume  fraction  and  polarizability  for  the 
ith  component,  respectively,  and  N  is  the  refractive  inaex. 

Ex-Situ  Ellipsometry  Studies 

If  the  InP  thermal  oxide  is  treated  as  one  homogeneous  film, 
i.e.  L1  =  o,  the  film  thickness (Lj)  and  refractive  index (N2)  can  be 
determined  from  the  measured  ¥  and  A  values  of  the  ex-situ 
experiment.  The  calculated  results  versus  the  etching  time  have 
been  plotted  in  Fig.  3  using  the  previously  obtained  value  3.521  + 
iO . 300  for  the  complex  refractive  index  of  InP  substrate (23) .  It 
is  clear  that  the  etch  rate  dramatically  changes  at  around  90 
seconds  with  the  averaged  etch  rates  before  and  after  that  point  of 
approximately  17  and  0.5nm/min,  respectively.  The  more  than  30 
times  change  in  the  etching  rate  indicates  a  significant  difference 
in  the  chemical  composition,  i.e.  a  two-layer  structure  for  the 
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oxide.  The  two-layer  oxide  structure  is  consistent  with  reported 
XPS  analyses (8 , 9 , 24 )  .  Although  the  exact  chemical  composition  of 
the  InP  thermal  oxide  is  still  unknown,  it  is  generally  accepted 
that  tl  a  oxide  is  composed  oi  an  outer  layer  of  a  mixture  of  ln203 
and  InP04  and  an  inner  layer  of  predominantly  InP04.  It  is  apparent 
that  the  composition  difference  causes  different  etch  rates  for  the 
InP  thermal  oxide  by  the  dilute  aqueous  HF  solution. 

Due  to  the  two-layer  structure,  the  one -film  model  used  for 
the  calculations  in  Fig.  3  is  valid  only  beyond  90  seconds.  From 
Fig.  3,  the  thickness  of  the  inner  la^erCI^),  of  12.2nm  car  be 
directly  obtained  as  the  thickness  at  the  turning  point,  while  the 
refractive  index  for  chis  layer  has  been  determined  to  be  2.09  by 
a/eraging  the  index  values  beyond  that  point  and  assuming  a  non¬ 
absorbing  oxide  for  this  wavelength,  i.e.  zero  imaginary  part  for 
the  index  (9)  .  Knowing  Lj  and  N2,  the  thickness  (L,y  and  refractive 
indexi’N,)  for  the  outer  layer  have  been  determined  to  be  I5.7nm  and 
1.94,  respectively,  according  to  equation  (1).  The  final  picture 
tor  the  thermal  oxide  is  given  in  the  insert  of  Fig.  3.  As 
mentioned  above,  the  one  film  model  is  only  valid  for  the  etch  time 
greater  than  90  seconds.  For  the  initial  region,  calculations 
based  on  the  two-film  model  have  been  carried  out.  Since  results 
for  the  thicknesses  of  these  points  are  only  slightly  different 
from  the  values  shown  in  the  Fig.  3  due  to  small  difference  in 
index,  we  will  noc  address  the  thickness  correction.  Bergignat  et. 
al.  have  measured  the  dielectric  functions  for  InP  thermal  oxide 
grown  in  the  250  -  550°C  temperature  range  using  spectroscopic 
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ellipsometry (9)  .  An  one  film  model  was  used  in  their  data  analysis 
for  the  thermal  oxides  grown  at  all  temperatures  except  550°C,  for 
which  they  used  a  two-film  model.  For  the  oxide  grown  at  4  50°C,  at 
the  wavelength  of  632. 8nm,  the  refractive  index  was  reported  to  be 
1.86,  which  is  lower  than  our  values  for  either  the  inner  or  outer 
layer  of  the  oxide  grown  at  440°C.  The  difference  may  arise  from 
the  10°C  difference  in  oxidation  temperatures  and/or  a  different 
value  used  in  the  data  analysis  for  the  complex  refractive  index  of 
InP  substrates.  However,  the  index  value  used  in  their  work  can 
not  be  explicitly  obtained. 

In-Situ  Ellipsometry  Studies 

The  experimental  ^  and  A  values  for  the  in-situ  ellipsometric 
measurement  of  the  thermally  grown  InP  oxide  etched  by  the  dilute 
HF  solution  are  shown  in  Fig.  4  along  with  calculated  values.  The 
calculated  V  and  A  values  assume  constant  refractive  indices  for 
films  and  substrate,  and  first  a  decreasing  film  thickness  of  the 
outer  layer  and  then  the  inner  layer  and  with  both  layers  having  a 
14nm  starting  thickness.  Since  the  etchant  contains  only  0.1 
percent  HF,  the  refractive  index  (1.33  25)  of  pure  H20  at  the 
wavelength  of  632. 8nm  has  been  used  for  the  dilute  HF  solution(25)  . 
The  experimental  starting  point  was  at  ¥  =  15.04°  and  A  =  68.32°, 
and  represented  the  sample  under  d.  i.  water  prior  to  the  HF 
etching,  with  a  14.4nm  inner  layer  and  12nm  outer  layer.  Another 
oxide  sample  grown  similarly  has  been  examined  using  the  same 
etching  procedure  with  a  2  minute  time  interval  between  successive 
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measurements  and  has  yielded  nearly  identical  results.  As  evident 
in  the  Figure,  after  the  start  point  the  measured  data  neither 
follow  the  calculated  curve,  nor  approach  the  point  for  a  bare  InP 
substrate.  It  is  apparent  that  a  different  mechanism  other  than 
simply  film  etching  is  involved.  However,  the  above  reported  ex- 
situ  results  show  the  expected  behavior  for  two-layer  film  etching. 
The  oxide  film  roughening  was  expected  to  some  degree  during  the 
prolonged  etching.  The  effects  of  a  rough  film  were  estimated 
assuming  a  rough  film  with  varying  void  percentage.  It  was  found 
that  the  extremely  rough  film  could  not  account  for  10  percent  of 
the  digression  of  the  measured  A  from  the  expected  A.  Thus,  during 
the  etching  the  sequence  of  events  must  include  the  formation  of 
layers  that  cause  large  changes  in  polarized  light  in  addition  to 
the  etching  of  the  two-layer  oxide  film.  In  a  separate  study,  a 
similar  behavior  was  observed  during  HF  etching  of  thermally  grown 
Si02  on  Si,  but  only  after  the  oxide  was  etched  away. 

Considering  the  results  above,  a  simple  model  for  this  new 
phenomenon  would  be  a  liquid  layer  sandwiched  by  the  residual  oxide 
and  the  bulk  etchant.  It  is  sensible  to  propose  that  there  is  a 
concentration  gradient  of  etch  products  in  the  liquid  layer 
extending  from  the  etching  oxide  layer.  With  these  assumptions, 
the  dielectric  response  of  this  liquid  layer  is  then  calculated 
based  on  the  Lorentz-Lorenz  relation (22)  which  is  appropriate  for 
individual  polarizable  species.  We  shall  proceed  with  the  details 
of  this  calculation. 

It  should  be  understood  that  single  wavelength  ellipsometry 
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has  restrictions  in  the  modeling  of  an  optical  system (26) .  First, 
the  dielectric  response  of  the  system  may  not  be  maximized  at  the 
wavelength  632.8  nm,  hence  the  sensitivity  may  not  be  optimum. 
Secondly,  by  fitting  a  single  measured  point  rather  than  a  set  of 
points,  as  would  be  the  case  for  spectroscopic  ellipsometry ,  the 
measurement  errors  can  propagate  through  the  modeling  to  the  final 
result.  As  a  consequence,  the  modeling  can  only  be  considered  to 
have  qualitative  significance.  However,  knowledge  gained  about  the 
sample  before  and  after  the  in-situ  etching  has  been  used  in  the 
modeling,  and  we  show  that  much  of  the  ambiguity  associated  wit 
single  wavelength  ellipsometry  can  be  eliminated. 

From  the  ex-situ  etching  experiment  above,  we  know  that  the 
etching  of  the  outer  layer  of  the  oxide  proceeded  relatively  fast. 
For  the  particular  sample  used  in  the  in-situ  measurement,  the 
thickness  of  the  outer  oxide  layer  was  about  12  nm.  This  thickness 
should  be  totally  removed  within  a  minute,  leaving  only  the  inner 
oxide  after  5  minutes  in  the  etching  environment.  The  refractive 
index  of  2.09  for  the  inner  oxide,  as  determined  in  the  above  ex- 
situ  measurement,  has  been  used  for  the  residual  inner  oxide.  In 
addition,  from  measurements  in  air  after  the  in-situ  experiment, 
the  thickness  for  the  final  residual  oxide  has  been  determined  to 
be  7.1  nm.  Using  the  two  film  model  and  starting  from  the  last 
measured  point  of  the  in-situ  data,  with  2.09  and  7.1  nm  as  the 
refractive  index  and  the  thickness  of  the  residual  oxide  on  InP  in 
the  dilute  HF  solution,  a  0.8  nm  unknown  top  layer  with  an  index  of 
3.07  is  calculated  to  be  between  the  bulk  solution  and  the  inner 


12 


oxide.  We  believe  that  this  layer  is  a  liquid  layer  because  it 
disappears  upon  removal  of  the  sample  from  the  solution  cell  and 
subsequent  d.i.  rinse  and  blow  dry.  From  this  last  point  towards 
the  beginning  of  the  etching  experiment,  we  assumed  that  the  inner 
oxide  is  thicker.  However,  any  significant  increase  in  the 
thickness (Lj)  beyond  7.1nm  did  not  result  in  an  acceptable  model 
until  the  etching  time  was  less  than  15  minutes.  Therefore,  the 
thickness  of  the  inner  oxide  has  been  fixed  at  Lj  =  7.1nm  in  the 
modeling  for  all  the  data  points  except  for  the  initial  two  where 
Lj  has  been  determined  to  be  8.1  and  7.4nm,  respectively,  for  the 
first  and  second  points  as  the  fitting  parameters  for  the  two  film 
model.  The  modeling  results  for  the  liquid  layer,  L1  and  N,,  are 
summarized  in  Fig.  5.  The  initial  5  minute  etch  resulted  in  a  121 
nm  liquid  layer  with  an  index  of  1.34,  and  the  thickness  of  the 
liquid  layer  decreased  quickly  to  about  1  nm  at  around  30  minutes 
while  the  index  for  the  layer  did  not  change  significantly.  Then, 
the  thickness  fluctuated  slightly  around  1  nm,  and  the  index  rose 
sharply  for  the  next  10  minutes  leveling  off  to  a  value  greater 
than  3.0  after  120  minute  etch.  Figure  6  shows  a  diagram  of  this 
etching  system,  the  liquid  layer  L1  and  oxide  layer  Lj. 

In  order  to  understand  the  nature  of  the  liquid  layer,  the 
chemical  reactions  involved  in  the  etching  process  should  be 
considered.  Since  the  HF  solution  can  etch  fused  silica  as  well, 
the  effects  of  the  reaction  of  the  HF  solution  with  the  cell  needs 
to  be  evaluated.  It  has  been  reported  that  the  etching  rate  of  the 
fused  silica  by  1:500  HF(49%) iHjO  was  about  0.3  nm/min(16) .  For 
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the  particular  cell (150cm3  in  capacity)  used  in  our  experiments, 
the  total  number  of  Si02  molecules  etched  away  from  the  cell  was 
calculated  to  be  1.3  X  1015  per  minute.  Assuming  that  the  etch 
products  were  SiF62‘  ions,  only  7.8  x  1015  F*  ions  were  consumed  in 
one  minute.  For  up  to  1000  minutes,  which  are  much  longer  than  any 
of  our  in-situ  measurements,  7.8  x  1018  F'  ions  were  consumed  and 
1.3  x  1018  SiF62*  ions  were  formed.  Compared  to  the  total  amount  of 
F*  ions  in  the  cell,  5.1  x  1021,  the  number  of  F'  ions  consumed  was 
negligibly  small.  The  SiF62*  ions  formed  accounted  for  only  0.3 
ppm.  Hence,  we  conclude  that  neither  the  concentration  nor  the 
physical  properties,  i.e.  refractive  index,  of  the  etching  solution 
should  be  significantly  affected  by  this  process. 

The  behavior  of  the  liquid  layer  as  shown  in  Fig.  5,  for  the 
in-situ  measurement  before  and  after  30  minutes  is  distinctly 
different.  Thus,  we  separately  examine  the  two  regimes.  To  start 
we  consider  that  for  the  InP  oxide  etched  by  the  aqueous  HF  three 
steps  were  involved:  1)  F*  ions  diffuse  to  the  oxide  surface,  2) 
F'  ions  react  with  the  outer  layer  of  the  oxide,  and  3)  the 
dissociated  species  diffuse  away  from  the  interface.  Initially, 
the  surface  reaction  proceeds  at  a  relatively  high  rate  as 
previously  discussed,  with  the  formation  of  significant  amounts  of 
ions  and  particles  at  the  oxide  surface  and  diffusing  into  the  bulk 
solution.  As  a  result,  a  concentration  gradient  is  formed  near  the 
solution  solid  interface,  which  in  turn  causes  the  optical  response 
in  the  in-situ  ellipsometric  measurement.  Some  possible  reactions 
taking  place  at  the  oxide  surface  are: 
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In203  +  6H*  =  2In3+  +  3H20 
InP04  +  2H+  =  In3+  +  H2P04' 

P205  +  3H20  =  2H2P04*  +  H* 

In3+  +  6F‘  =  InF63' 

A  mixture  of  the  ions  such  as  IT",  F*,  In3-*,  H2P04',  InF63‘,  and  solid 
particles  such  as  InP04  and  P205  in  the  vicinity  of  the  oxide,  can 
be  imagined  to  constitute  a  concentration  gradient  in  L1  extending 
outward  from  oxide  surface  depicted  in  Fig.  6.  If  all  the  P  and  In 
species  in  the  etch  products  are  counted  in  the  form  of  H3P04  and 
InF3,  there  would  be  approximately  1.42  and  4.02  x  1016  molecules, 
respectively,  formed  after  the  initial  5  minute  etch.  Assuming 
that  these  species  are  uniformly  distributed  within  a  distance  from 
the  solid  surface,  and  beyond  the  distance,  the  concentration  of 
the  diffused  species  becomes  negligible,  an  abrupt  boundary  can  be 
defined  as  a  plane  parallel  to  the  solid  surface  with  a  vertical 
distance  as  the  measured  thickness  for  the  liquid  layer  and  an  area 
of  the  sample  size,  i.e.  1  cm2.  Using  the  Lorentz -Lorenz  relation, 
the  volume  fractions  for  H3P04  and  InF3  have  been  calculated,  where 
5.3  x  10*24  cm3  has  been  used  for  the  polarizability  of  H3P04(27). 
We  found  no  reference  for  the  polarizability  value  of  InF3,  so  it 
was  treated  as  a  parameter.  For  these  specific  calculations,  4.4 
x  10*24  cm3  has  been  used,  and  it  was  found  that  the  numerical  value 
did  not  qualitatively  affect  the  calculations.  Relative  to  the 
total  amount  of  the  P  and  In  species  formed,  the  percentage  of  the 
moleculars  remaining  in  the  boundary  has  also  been  calculated. 
Similar  treatments  have  been  carried  out  for  several  other  points 
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and  summarized  in  Table  1.  The  first  row  of  the  table  lists  the 
etching  time,  and  the  second  and  the  third  are  the  thicknesses  and 
refractive  indices,  respectively,  for  L1  from  Fig.  5.  The  total 
amounts  of  P  and  In  species  produced  by  the  etch  in  the  forms  of 
H3P04  and  InF3  are  listed  in  the  fourth  and  fifth,  respectively. 
These  values  were  obtained  by  counting  the  numbers  of  the  ln203  and 
InP04  molecules  in  the  oxide  before  they  were  etched  away.  The 
bottom  row  is  the  percentage  of  total  etch  products  remaining  in 
the  liquid  layer.  As  shown  in  the  table,  the  refractive  indices 
and  volume  fractions  for  both  species,  did  not  change  much  before 
25  minutes  of  etching  time,  while  the  thickness  and  the  percentage 
of  etch  products  remaining  in  the  layer  decreased  steadily.  The 
percentage  of  the  remaining  etch  species  and  liquid  layer  thickness 
from  Table  1  were  plotted  versus  the  etching  time  in  Fig.  7a  and 
7b,  respectively.  Fig.  7a  shows  a  linear  decrease  of  the  total 
etch  products  remaining  in  the  liquid  layer,  while  Fig.  7b  shows 
that  the  liquid  layer  decreases  linearly  from  10  to  20  minutes, 
then  falls  at  an  even  faster  rate.  The  Maxwell -Garnett  effective 
medium  approximation  treatment,  in  which  the  etch  solution  was 
considered  as  a  host  medium (H20)  containing  dissociated  particles 
(H3P04  and  InF3)  ,  has  also  been  carried  out  to  examine  the  liquid 
layer  and  similar  to  the  Lorentz -Lorenz  results  we  found  a  linear 
decrease  in  the  percentage  of  both  total  P  and  In  remaining  in  the 
liquid  layer.  However,  the  value  for  the  first  point  and  the  slope 
of  the  linear  decrease  were  much  lower  than  that  from  Lorentz- 
Lorenz  treatment.  With  this  modeling,  the  physical  picture  of  this 
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particular  etching  phenomenon  becomes  clear.  Initially,  the  surface 
reaction  of  F'  with  the  oxide  proceeds  relatively  fast,  and  the 
etch  products  diffusing  away  from  the  oxide  surface  form  a  thick 
liquid  layer  in  a  short  time.  After  the  removal  of  the  outer  oxide 
layer,  the  reaction  slows  down  and  eventually  stops,  and  so  does 
the  supply  of  the  P  and  In  species.  Thus,  the  out  diffusion  of 
these  species  becomes  dominant  due  to  the  circulation  of  the  HF 
solution,  which  maintains  a  constant  concentration  in  the  bulk 
solution,  and  transports  the  out-diffused  species  out  of  the  cell. 
The  liquid  layer  starts  to  shrink  slowly  and  then  at  a  linear  rate. 
After  20  minutes,  as  shown  in  Fig.  7b,  the  liquid  layer  thickness 
drops  at  a  higher  rate.  The  linear  decrease  of  the  thickness  from 
10  to  20  minutes  can  be  described  by: 

L  -  169  -  5t ;  (4) 
where  L  is  the  thickness  of  the  liquid  layer  in  nanometers,  and  t 
is  the  etching  time  in  minutes.  Equation  (4)  also  represents  the 
rate  at  which  the  etch  products  diffuse  out  of  the  surface  region. 

So  far,  we  have  modeled  the  in-situ  ellipsometric  data  for  HF 
etching  InP  oxide  and  established  a  simple  physical  picture  for  the 
complex  etching  system.  If  we  keep  a  constant  refractive  index  for 
the  liquid  layer,  1.34  in  this  case,  vary  the  thicknesses  according 
to  equation  (4) ,  and  use  the  same  inner  oxide  layer  as  described 
above  and  the  two  film  model,  the  ¥  and  A  can  be  obtained,  in  Fig. 
8,  curve  b.  The  calculated  data  for  a  two  layer  oxide  in  H20  of 
Fig.  4  is  shown  again  in  Fig.  8  as  curve  a.  Before  the  introduction 
of  HF  solution,  the  measured  ¥andA  point  was  close  to  curve  a,  then 
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the  measured  data  jumped  onto  curve  b(see  Fig.  4)  and  followed  it 
very  closely  up  to  30  minutes  in  the  etching  solution.  Afterwards, 
the  measured  data  departed  from  curve  b.  The  deviation  can  be 
accounted  for  by  considering  the  contraction  of  the  liquid  layer 
and  will  be  discussed  below. 

Before  we  proceed  to  the  next  regime  of  the  in-situ 
measurement,  we  need  to  examine  and  evaluate  the  assumptions  and 
approximations  made  in  the  above  discussion.  Firstly,  we  have 
assumed  that  an  abrupt  boundary  layer  existed.  Within  the  boundary 
layer,  the  ions  and  solid  particles  were  uniformly  distributed, 
while  outside  it,  there  were  negligible  amounts  of  ions  and 
particles.  Therefore,  an  isotropic  and  homogeneous  optical  layer 
could  be  assumed  in  the  two  film  model.  This  is  obviously  only  a 
rough  approximation  of  a  continuous  concentration  gradient. 
However,  because  of  the  low  concentrations  and  the  gradual  change 
from  the  oxide  surface  to  the  bulk  solution,  the  optical 
properties,  i.e.  refractive  index,  for  both  sides  of  the  boundary 
were  only  slightly  different,  we  believe  that  this  is  a  reasonable 
approximation  for  this  particular  system.  Secondly,  the  P  and  in 
species  formed  were  accounted  for  as  only  H3P04  and  InF3,  with  the 
estimated  polarizability  values  for  H3P04  and  InF3  having  been  used 
in  the  calculation.  In  reality  a  large  number  of  species  may  exist 
in  solution  and  each  individual  species  is  expected  to  be  different 
in  optical  respones.  Nevertheless,  the  physics  and  the  qualitative 
meanings  of  our  simple  treatment  is  a  useful  average  way  to  treat 
the  complex  system.  The  third  assumption  was  that  we  neglected  the 


surface  roughness,  surface  electrical  charge  distribution,  ions 
adsorbed  by  the  oxide,  as  well  as  alignment  of  the  H20  molecules 
along  their  permanent  dipole  moments  in  the  vicinity  of  the  solid 
surface.  All  these  factors  could  affect  the  observed  results  in 
one  way  or  another,  but  apparently  from  the  success  of  the  simple 
model  none  of  them  predominated. 

Beyond  30  minutes,  a  different  mechanism  was  dominant.  As  the 
thickness  of  the  liquid  layer  decreased  quickly  to  about  1  nm,  the 
refractive  index  markedly  increased  and  then  leveled  off  as  shown 
in  Fig.  5.  The  rapid  contraction  of  the  liquid  layer  suggests  that 
certain  ions  and/or  molecules  have  been  driven  towards  the  solid 
surface,  and  formed  a  dense  liquid  layer.  The  layer  condensed 
slowly  for  about  95  minutes  and  eventually  saturated.  Curve  c  of 
Fig.  8  are  calculated  data  for  a  1  nm  liquid  layer  on  the  residual 
oxide  with  a  varying  index  for  the  liquid  layer.  The  fit  of  this 
curve  to  the  rest  of  the  measured  data  supports  our  ellipsometric 
modeling  for  this  regime,  as  seen  by  the  experimental  data  in  Fig. 
4  with  the  model  in  Fig.  8. 

The  driving  force  for  the  condensation  is  likely  electrical  in 
nature,  and  possibly  results  from  either  the  net  electrical  charges 
in  the  residual  oxide  due  to  the  removal  of  the  counter  charges 
resident  in  the  outer  oxide,  or  the  specific  ion  adsorption  by  the 
inner  oxide,  or  both.  Applying  the  previous  results  from  the  two 
film  modeling  to  the  electric  double  layer,  the  thickness  and 
refractive  index  are  found  to  be  1.0  nm  and  3.1,  respectively.  As 
for  comparision,  the  thickness  and  refractive  index  of  the  double 
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layer  in  the  second  measurement  are  determined  to  be  1.2nm  and  2.8, 
respectively.  Stedman(12)  has  reported  an  investigation  of  the 
electric  double  layer  formed  at  a  metal  electrolyte  interface  using 
ellipsometric  techniques,  and  a  simple  theoretical  treatment  was 
also  discussed.  The  optical  response  of  the  components  of  the 
electric  double  layer  was  examined  separately  in  that  work.  The 
inner  layer  consisted  of  compressed  solvent  dipoles  that  resulted 
in  a  larger  refractive  index  than  the  bulk  solvent.  The  ions 
adsorbed  on  the  solid  surface  and  the  ions  in  the  Gouy  Chapman 
diffuse  layer  also  contributed  to  the  observed  optical  response. 
The  type  of  ions  and  their  distribution  in  the  electric  double 
layer  determined  the  dielectric  properties.  In  principle,  the 
refractive  index  can  be  estimated  using  the  Lorentz -Lorenz  equation 
if  the  concentration  of  the  adsorbed  ions  and  ion  distribution  in 
the  diffuse  layer  is  known.  As  discussed  above,  the  system  of  the 
HF  etching  InP  oxide  is  complicated.  A  variety  of  ions  were  formed 
and  contributed  to  the  electric  double  layer.  However,  the  detail 
knowledge  about  the  structure  of  the  electric  double  layer,  which 
is  essential  for  the  estimation  of  the  refractive  index,  could  not 
be  obtained  by  the  ellipsometric  technique  alone.  It  would  be 
valuable  to  combine  electrochemical  techniques  for  the  study  of  the 
solution-solid  interfacial  phenomena. 

Conclusions 

The  InP  oxide  surface  reaction  with  a  dilute  HF  solution  has 
been  studied  by  both  ex-situ  and  in-situ  ellipsometric  techniques. 
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The  two-layer  structure  for  the  thermally  grown  InP  oxide  has  been 
clearly  identified  and  characterized.  The  outer  layer  of  the  oxide 
which  has  a  lower  refractive  index  was  etched  much  faster  than  the 
inner  layer  which  has  a  higher  index  value. 

A  two-layer  optical  model  has  been  used  to  analyze  the  in-situ 
etching  ellipsometry  data.  A  liquid  layer  consisting  of  the  etch 
products  of  P  and  In  species  has  been  detected  for  the  initial 
regime.  Based  on  the  Lorentz-Lorenz  relation,  a  microscopic  model 
has  been  established  to  describe  the  liquid  layer.  A  linear 
process  can  account  for  the  out  diffusion  of  the  etch  species, 
which  forms  the  liquid  layer.  The  theory  agrees  qualitatively  with 
our  experimental  results.  An  electric  double  layer  was  observed  to 
form'  at  the  oxide-solution  interface.  The  thickness  and  refractive 
index  for  the  double  layer  has  been  estimated  as  1.0  nm  and  3.1, 
respectively.  With  the  sensitivity  of  ellipsometry,  in-situ 
solution  ellipsometry  is  a  promising  technique  for  the  study  of  the 
solid-solution  interface  phenomena. 
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List  of  Figures 


Figure  1:  A  sketch  of  the  automatic  ellipsometer . 

Figure  2:  The  illustration  of  the  fused  silica  solution  cell  for 
the  in-situ  ellipsometric  measurement.  The  ball  socket  joints  are 
required  in  order  to  align  the  entrance  and  exit  windows. 

Figure  3:  The  thicknesses  and  refractive  indices  for  the  thermal 
InP  oxide,  etched  by  the  dilute  HF  solution,  as  measured  by  the  ex- 
situ  ellipsometry. 

The  insert  shows  the  two  layer  structure  for  the  oxide. 

Figure  4:  The  measured  and  simulated  data  for  the  in-situ  etching 
in  'J'andA  space.  For  the  measured  data,  the  time  interval  between 
each  consecutive  points  is  5  minutes,  while  for  the  simulated  data, 
the  distance  between  each  two  points  represents  1  nm  increment  in 
film  thickness. 

Figure  5:  The  thicknesses  and  refractive  indices  for  the  liquid 
layer  resulted  from  the  two  film  modeling  of  the  in-situ  etching 
measurement  for  the  InP  oxide. 

Figure  6:  A  representation  of  the  liquid  layer  L.,,  which  consists 
of  a  variety  of  ions  and  neutral  particles,  formed  on  top  of  the 
thermal  InP  oxide,  1^,  during  the  etching  process. 

Figure  7:  a)  the  percentage  remaining  of  the  etch  products  in  the 
surface  region  versus  the  etching  time  from  the  calculation  of  the 
Lorentz-Lorenz  relation,  and  b)  the  liquid  layer  thickness  resulted 
from  the  two  film  modeling  versus  the  etching  time. 
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Figure  8:  The  in-situ  measured  data (open  circles)  and  three 
simulated  curves  labelled  a,  b,  and  c.  Curve  a  is  calculated 
assuming  a  two  layer  oxide  as  in  Figure  4 ,  curve  b  assuming  a 
constant  index  and  a  linearly  changed  thickness  for  the  liquid 
layer,  and  curve  c  for  a  1  nm  liquid  layer  with  variable  index. 

List  of  Tables 

Table  1:  A  summary  of  the  calculations  of  the  volume  fractions  for 
P  and  In  species  in  the  liquid  layer  for  the  etching  time  5  to  30 
minutes.  All  the  P  species  have  been  counted  in  the  form  of  H3P04, 
and  In  species  in  InF3. 
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Table  1 


Time  (Min) 

5 

10 

15 

20 

25 

30 

Exp .  Thickness 
(nm) ,  L, 

121 

119 

96 

69 

19 

1.2 

1.34 

1.34 

1.34 

1.34 

1.36 

2.27 

1.42 

1.49 

<N 

in 

• 

H 

1.52 

1.52 

1.52 

InF,  (xlO16) 

4.02 

4.14 

4.19 

4.19 

4.19 

4.19 

Remaining  (%) 

98 

90 

71 

51 

20 

13 
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